log_regularized_beta Function

public pure elemental function log_regularized_beta(a, b, x, max_iterations) result(rst)

Log regularized beta; invalid inputs or nonconvergence return NaN. The leading factor and complementary branch are evaluated in log space.

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(in) :: a

Finite, positive beta shape parameters.

real(kind=real64), intent(in) :: b

Finite, positive beta shape parameters.

real(kind=real64), intent(in) :: x

Integration limit on [0, 1].

integer(kind=int32), intent(in), optional :: max_iterations

Continued-fraction iteration limit; default 1000.

Return Value real(kind=real64)

Log ratio; x = 0 gives -infinity, x = 1 gives zero.